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Using a real-space renormalization-group approximation, we study the anisotropic quantum 
Heisenberg model on hierarchical lattices, with interactions following aperiodic sequences. Three 
different sequences are considered, with relevant and irrelevant fluctuations, according to the Luck- 
Harris criterion. The phase diagram is discussed as a function of the anisotropy parameter A (such 
that A = and A = f correspond to the isotropic Heisenberg and Ising models, respectively). We 
find three different types of phase diagrams, with general characteristics: the isotropic Heisenberg 
plane is always an invariant one (as expected by symmetry arguments) and the critical behavior 
of the anisotropic Heisenberg model is governed by fixed points on the Ising-model plane. Our 
results for the isotropic Heisenberg model show that the relevance or irrelevance of aperiodic mod- 
els, when compared to their uniform counterpart, is as predicted by the Harris-Luck criterion. A 
low-temperature renormalization-group procedure was applied to the classical isotropic Heisenberg 
model in two-dimensional hierarchical lattices: the relevance criterion is obtained, again in accor- 
dance with the Harris-Luck criterion. 

PACS numbers: 



I. INTRODUCTION 

The investigation of systems displaying inhomogeneous 
or disordered interactions is an active field of research 
[? ? ]. From the experimental point of view, many 
of the materials found in nature come with impurities; 
also, modern techniques are able to build materials with 
controlled composition, such that two or more different 
atoms are combined in a given order. Theoretically, one 
may be concerned with possible changes on the critical 
behavior of systems, with the introduction of random dis- 
order or inhomogeneuos deterministic interactions, when 
compared to their homogeneous counterpart [? ]. For 
quenched random disorder, the Harris criterion [? | 
states that, if the pure-system's specific-heat exponent, 
a, is positive (negative), the critical behavior of the disor- 
dered model is diffferent from (the same as for) the pure 
model. 

The discovery of quasi-crystals [? ] has motivated an 
intense research on the behavior of models with interac- 
tions following aperiodic sequences [? ]: numerical [? ] 
as well as analytical results [? ? ] have been obtained. 
Many works concentrated on classical models, like Ising 
and Potts [? ? ] ones, but some attention has been drawn 
to quantum models in one dimension [??????]. In 
Ref. ? the ground-state properties of a two-dimensional 
quantum model have been analyzed. However, no work 
has focused on the finite-temperature critical behavior of 
quantum models in dimensions two or above, to the best 
of our knowledge. 
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A convenient model to address the role played by quan- 
tum effects and aperiodicity on critical phenomena is 
the anisotropic Heisenberg one, with interactions follow- 
ing aperiodic and deterministic sequences. On the other 
hand, the way these sequences are constructed and the 
idea behind renormalization-group calculations make hi- 
erarchical lattices a natural choice for the study. There- 
fore, in this work we treat the anisotropic ferromagnetic 
Heisenberg model on four different hierarchical lattices, 
with different Hausdorff fractal dimenions. Three dif- 
ferent aperiodic sequences are treated, corresponding to 
bounded and unbounded fluctuations. We will study here 
mainly the relevance of the introduction of aperiodicity, 
in the renormalization-group sense. This aspect is gener- 
ally addressed by the so-called Harris-Luck criterion [? ] . 
According to this criterion, the relevance of a given ape- 
riodic sequence is connected to the crossover exponent, 
4>, given by [? ? ] : 

<j> = l-d a v Q (i-u), (l) 

where d a is the dimension the aperiodic sequence acts on, 
Uq is the correlation lenght's critical exponent of the pure 
model and uj is the fluctuation exponent. This exponent 
is defined through g ~ N u , where g is the fluctuation in 
the number of a given letter of the sequence (below, we 
discuss this point further). For cj> > 0, the critical behav- 
ior of the aperiodic and uniform models are in different 
universality classes. For <fr < 0, both aperiodic and uni- 
form models have the same set of critical exponents. For 
the marginal case, = 0, critical exponents depend on 
the ratio between the two interaction constants [? ? ]; 
we will not duscuss this case further. Our choices of hier- 
archical lattices and aperiodic sequences allow for values 
of (f> greater os smaller than zero, as well as for differ- 
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ent fractal dimensions of the lattice. The aperiodicity is 
chosen such that d a = 1 in all cases we study. 

This work is organized as follows. In the next section 
we review some basic concepts of aperiodic sequences, 
which will be important to our work. In Scction HTTl we de- 
fine the model and outline the real-space renormalization- 
group approach we use and in Section IIVI we present our 
results. In Section [V] we discuss and summarize the re- 
sults. 



II. APERIODIC SEQUENCES 

The aperiodic sequences used in this work are obtained 
by the iteration of substitutuon rules working on an al- 
phabet. Each letter of the alphabet is replaced by a se- 
quence of letters and words are formed at each stage. We 
will be concerned with a two- letter alphabet, A and B, 
and the usual convention is that the initial word is A. 
More specifically, we will study the following aperiodic 
sequences: 

(i) A — > AB; B — > A A, i.e., from a given word of the 
sequence, the next word is obtained by substituting 
AB for every letter A and AA for every letter B in 
the previous word. The first stagess of this sequence 
are A — > AB — ► ABA A -> ABAAABAB 
This is the so-called period-doubling sequence. 



(ii) A -> ABB; B - 
sequence are: A 



AAA. The first stages of this 
-> ABB -> ABBAAAAAA -> 



(iii) A -> AAB; B -> AAA. In this case, the 
first words of the sequence are: A —> AAB — > 
AAB AAB AAA 

The geometrical characteristics of these sequences are 
obtained from the substitution matrix Ai, which con- 
nects the number of letters A and B after one application 
of the iteration rule, namely: 



N {n + l) 



= M 




(2) 



Therefore, after n iterations the total number of letters 
in the word, N^ n \ is given by: AT(«) = N ( A n) + N ( ^ ] = 
M n N^\ where N A n) and are the number of letters 
A and B after n iterations, respectively, and A^ ' = 1 

and N B °^ — for the initial word. 

The substitution matrices for the three aperiodic se- 
quences defined above are: 

(i)A4=(j o 



(ii) Ai 



1 3 

2 



The total number of letters grows exponentially with 
the number of iterations n: 



N 



A™, n — ► co, (3) 

where N = lim„^oo N^ n ' and Ai is the greater eigenvalue 
of Ai. For the three sequences studied in this work, this 
exponential growth is valid for all n and Ai = 2 for se- 
quence (i) and Ai = 3 for sequences (ii) and (iii). One 
can define the fluctuation in a given letter, say A, as 
(™) — TV 
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(») 



PaN( u \ where pa is the fraction of let- 
ters A in the infinite word, i.e, after n applications of the 
iteration rules, with n — > co. The fractions pa and pb 
are proportional to the first and second entries, respec- 
tively, of the eigenvector corresponding to the greater 
eigenvalue. It is possible to show that: 



(4) 



where g = limn^oo and A2 is the smaller eigenvalue 
of Ai. Therefore, using Eqs. ([3]) and Q, one can show 
that: 
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u> 



In 1 A3 1 
I11A1 



(5) 



The exponent u> is crucial for the crossover exponent, 
as outlined in the previous section. Its value is w = 0, 
ln(2)/ln(3) and for sequences (i), (ii) and (iii), respec- 
tively, as can be easily calculated from their substitution 
matrices. We will disccus the results for the crossover 
exponent in Section HVl 



III. MODEL AND FORMALISM 

The reduced Hamiltonian of the anisotropic Heisenberg 
model is given by: 

~m=Y, K v K 1 - A )«^ + + < ( 6 ) 

<i,3> 

where (3 = 1/ksT, fc^ being the Boltzmann constant and 
T the temperature, af is the component a of a spin-1/2 
Pauli matrix on site j,0<A<1 (A = corresponds to 
the isotropic Heisenberg model and A = 1 to the Ising 
model), the sum is over all first-neighbor bonds on a hi- 
erarchical lattice and the exchange constants Kij = (3Jij 
follow an aperiodic sequence in a given direction of the 
lattice. See Figs. []] and [2] for examples of hierarchical 
lattices with b — 2 and 6 = 3, respectively: these lattices 
are built of bonds connected in parallel, each one 
consisting of b bonds connected in series, where d is the 
fractal dimension of the lattice. In this work, we have 
treated lattices with d = 2 and d = 3. 

We use a real-space renormalization-group approach; 
a partial trace is taken over internal spins on suitable fi- 
nite lattices and a renormalizcd Hamiltonian is obtained, 
namely: 



exp {K[ 2 [(1 - A')(of of + a\al) + afa$\} 
Tr W e-^, (7) 
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(a) 



(b) 



FIG. 1: Hierarchical lattice with b — 2, suitable for the study 
of the period-doubling sequence (sequence (i); see text). We 
show the renormalization for the coupling constant Ka (part 
(a)) and for the coupling constant Kb (part (6)). Note that 
the construction of the hierarchical lattice is made in the re- 
verse order of the renormalization-group procedure. 





FIG. 2: Hierarchical lattice with 6 = 3, suitable for the study 
of sequences (ii) and (iii) (see text). We show the renor- 
malization for the coupling constant Ka (part (a)), where 
Ki — Kb for sequence (ii) and Ki = Ka for sequence (iii), 
and for the coupling constant Kb (part (b)). 



where K' 12 = K A in Figs. [lja) and[2la) and K' 12 = K B 
in Figs. QJ&) and[^&), Trr^ is a partial trace, taken over 
all spins in Figs. Q]and[2l except <j\ and 02, and (¥H is 
the reduced Hamiltonian of the cell on the left-hand sides 
of Fig. □ (a) and Q] (b) and Fig. (a) and EH (b). The 
method to calculate the partial trace in Eq. {7J for quan- 
tum systems was introduced in Ref. ? and an important 
simplification was proposed in Ref. ? , where the whole 
process is explained in detail. This method has been suc- 
cessfully applied in the study of ferromagnetic, antifcr- 
romagnetic and spin-glass quantum models. The formal- 
ism is specially suitable to obtain multidimensional phase 
diagrams and qualitative results, indicating universality 
classes and possible crossover phenomena. It is worth 
mentioning that, although A is initially uniform and the 
aperiodicity acts only on the interaction parameter K , 
after the first iteration of the renormalization-group the 
anisotropy is no longer the same for every bond. This 
fact has to be taken into account when deriving the re- 
cursion relations. We refer the reader to Refs. ? and ? 
for details. 



IV. RESULTS AND DISCUSSION 

Some features are common to all three sequences: the 
isotropic Heisenberg model (A = 0) is an invariant sub- 
space, in the renormalization-group sense. The same ap- 
plies for the Ising model (A = 1). The reason is that 
the introduction of exchange constants with different val- 
ues does not change the symmetry of these two models. 
Therefore, the flow will not leave the corresponding sub- 
spaces. Moreover, the critical behavior for < A < 1 
will be determined by the stability of the Ising-modcl 
non-trivial fixed points, i.e, the flow for any initial value 
of A in that range is towards the A = 1 subspace. The 
Ising model with aperiodic interactions has already been 
treated and our results for A = 1 agree with those in 
Ref. ? . Particularly, the stability of the uniform fixed 
points, with respect to the introduction of aperiodicity, is 
in accordance with the Harri-Luck criterion. Therefore, 
from now on we will restrict ourselves to the isotropic 
Heisenberg model subspace. 

In the renormalization-group framework, the stability 
of fixed points is given by the eigenvalues of the matrix of 
the linearized renormalization-group equations (LRGE) 
[? ]. Since the relevant fixed points (filled squares in 
Fig. [3J are the non-trivial ones, one of these eigenvalues 
is always greater than one and corresponds to the flux 
along the uniform sub-space (traced lines in Fig. [3]). In 
this figure we show the qualitative picture we obtain for 
irrelevant (part (a)) and relevant (part (b)) aperiodic se- 
quences for d = 3. The uniform model corresponds to 
the straight line at 45°, where Ka = Kb; the fixed point 
is always unstable along this line. The relevance of the 
aperiodicity is given by the stability along the other di- 
rection (continuous lines leaving the fixed points in Fig. 
[3J: in part (a) of the figure, the aperiodicity does not 
change the critical behavior, compared to the uniform 
model, while in part (b) a new universality class emerges 
when aperiodicity is introduced. Note that the phase di- 
agrams for lattices with d = 2 are qualitative different 
from the ones in Fig. [3] the non-trivial fixed points are 
at zero temperature (K = 00) and, therefore, the "aperi- 
odic" direction is not physically accessible in these cases. 
Nevertheless, the relevance of the aperiodicity is correctly 
described by the renormalization-group formalism, as we 
will see below. 

A technical point is worth mentioning here. For all 
sequences we treat in this work the structure of the ma- 
trix of the LRGE, evaluated at the uniform fixed point, 



K* = K* 



( 


dK A 




dK A 




dK A 










K* 










ok' b 


V 


dK A 


K* 


8Kb 



A"* 



A"* 



Since: 



OKg/dK, 



K> 



A* 



dK A /dK B 



A* 



(8) 
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the eigenvalues of the above matrix are: 



dK A 



A* 



dK f 



(9) 



A ' 



The former corresponds to the uniform model, as dis- 
cussed previously, and, therefore, is always greater than 
1. The absolute value of the latter eigenvalue determines 
the relevance of the aperiodicity. 



(b) 



FIG. 3: Qualitative phase diagram for hierarchical lattices 
with d — 3. (a) Renormalization-group flux for irrelevant 
aperiodic sequences; (b) Renormalization-group flux for rel- 
evant aperiodic sequences. The straight lines at 45 degrees 
represent the uniform model, Ka = Kb, and the fixed point 
is relevant along this direction. 

For classical models, like the Ising and Potts ones [? ], 
the matrix of the LRGE is proportional to the transpose 
of the substitution matrix (see Eq. ([2|) above), M. . This 
property still holds true for quantum models and for ape- 
riodic sequences such that the hierarchical lattices which 
renormalizes into K A is symmetric with respect to the ex- 
change of the iterations Ka and Kb, as in part (a) of Fig. 
[TJ This brings an important simplification for the calcula- 
tion of the former matrix, since the proportionality factor 
can be obtained from the recursion relation for the uni- 
form model, i.e, from the renormalization of Kb- Then, 
one can obtain the other two elements, dK A / 8Ka 



and dK A /dKs 



from the substitution matrix, 



K* 

with 



A'* 



no need to work out the recursion relation for the ape- 
riodic model (this relation involves bonds with different 
values, which may make it cumbersome to be calculated). 

Note that, while the proportinality between the matrix 
T 

of the LRGE and M. is always true for classical models, 
it fails for quantum models when the lattice is not sym- 
metric with respect to the exchange of the iterations Ka 
and Kb (like the one in part (a) of Fig. [21 for example) . 

For the hierarchical lattices with b = 3, we have to 
resort to the Migdal-Kadanoff approximation, since the 
number of sites involved in the renormalization is rather 
large. This approximation is equivalent to treating the 
cell "by pieces", i.e., to renormalizating first the b bond in 
series and then combining the renormalized bonds 
in parallel. An alternative (and more precise, since com- 
mutation aspects are taken into account at the cell level) 
procedure is to renormalize the whole cell. This point is 



extensively discussed in Ref. ? , where it is shown that 
treating the cell "by pieces" gives good qualitative results 
for ferromagnetic quantum models, when compared to 
renormalizing the whole cell. The comparison between 
these two procedures has also been made for antiferro- 
magnetic quantum models [? ] and again the agreement 
is quite good (note, however, that for models such that 
frustration effects are present, the two approaches give 
different qualitative results [? ])). In this work, we have 
applied the two procedures outlined above for 6 = 2; the 
qualitative (and sometimes even the quantitative) agree- 
ment is excellent and put the Migdal-Kadanoff approxi- 
mation made for b = 3 on a firmer basis. 

The hierarchical lettices we treat here are those with 
b = 2, d = 2 and 3, (see Fig. [JJ and b = 3, d = 2 and 3 
(see Fig. [5]) . The corresponding results are: 

(a) b = 2, d = 2: in this case, the results are the same 
treating the cell as a whole or within the Migdal- 
Kadanoff approximation. The critical temperature, 
T c , and the correlation-length's critical exponent, 
v, are known exactly for the two-dimensional uni- 
form model, namely T c = [? ] and v = oo. 
We obtain these exact results with our procedure. 
Therefore, according to the Harris-Luck criterion, 
Eq. flj), the crossover exponent is negative for any 
sequence. In fact, the eigenvalues for this model 
are: Ai = 1 and A2 = —1/2. The former corre- 
sponds to the pure-model's critical behavior; the 
second determines the irrelevance of the aperiodic- 
ity, as predicted by the Harris-Luck criterion. The 
negative sign of the second largest eigenvalue is a 
signature of aperiodic systems [? ]; 

(b) b = 2, d — 3: in this case, we have 7) renormal- 
ized the cell as a whole and 77) used the Migdal- 
Kadanoff approximation, as in the last item. In 
both procedures we obtain a finite critical temper- 
ature (T c = 2.70 for the former and T c = 2.91 for 
the latter), as expected. The approximated value 
for the critical exponent v is: 7) 1.511, and 77) 
1.398. Since w = for the period-doubling se- 
quence, the crossover exponent is negative in both 
approximations, namely: 7) <f> — —0.511, and 77) 
<p = —0.398. Therefore, we expect the aperiodicity 
to be irrelevant. The eigenvalues of the matrices of 
the LRGE are: 7) Ai = 1.582; A 2 = -0.791; and 
77) Ai = 1.642; A 2 = -0.821. As expected, in both 
cases the aperiodicity is irrelevant and the smaller 
eigenvalue is nagative, as commented above; 

(c) b = 3, d = 2: again, T c = and v = 00 for the uni- 
form model. Our procedure (see Fig. [5]) obtains the 
correct critical temperature but the value obtained 
for v is an excellent approximation but not the ex- 
act one, namely v — 8.494. As discussed above, 
we have to resort to the Migdal-Kadanoff proce- 
dure, in this case, but, in view of the comparison 
made in (o) and (&), we believe that the physical 
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behavior is correctly described by this approxima- 
tion. For this lattice, we have studied sequences (ii) 
and (Hi) (defined in Section [Til- For sequence (ii), 
the wandering exponent is u> = hi(2)/ ln(3) and the 
crossover exponent is 4> = —2.135, i.e., the aperi- 
odicity is irrelevant, according to the Harris-Luck 
criterion; the situation is analogous for sequence 
(Hi), where u> = and then <fi = —7.494: the ape- 
riodicity defined by this sequence is also irrelevant. 
For the eigenvalues of the substitutional matrices 
we obtain: (ii): A 1 = 0.879; A 2 = -0.654; (Hi): 
Aj = 0.879; A 2 = -0.224. Therefore, the aperi- 
odicity is irrelevant for both cases, as predicted by 
the Harris-Luck criterion. Note that, since the rel- 
evant fixed point is at K\ = K B = oo, the value 
Ai = 0.879 means that this point unstable, in the 
renormalization-group sense; 

(d) b = 3, d = 3: we have used the Migdal-Kadanoff ap- 
proximation and obtained T c — 1.92 and v — 1.551 
for the uniform case. Recalling the values for the 
wandering exponents for sequences (ii) and (Hi) 
(see (c) above), the crossover exponents are given 
by (j> = 0.427 and <j) = -0.551, respectively. The 
eigenvalues of the matrix of the LRGE are: (ii): 
Ai = 2.030; A 2 = -1.363; (Hi): A 1 = 2.030; 
A 2 = —0.667. Therefore, sequence (ii) is rele- 
vant and sequence (Hi) is irrelevant, again in ac- 
cordance with the Harris-Luck criterion. For the 
relevant fixed cycle of period two emerges, 

as already found in Ref. ? . The location of 
this stable two-cycle is (K\ = 0.417, K* B = 3.90); 
(K* A = 1.85, K* B = 0.325). The eigenvalues of 
the LRGE associated with this double iteration are 
Ai = 3.89; A 2 = 0.250. The specific heat critical 
exponent associated with the aperiodic fixed cycle, 
a a , is calculated from the relation: 



2-d- 



lnfe 2 



where d is the fractal dimension of the hierarchical 
lattice and b is the scaling parameter associated 
with onde iteration. In this case, we obtain a a = 
—2.855, which is smaller than its counterpart for 
the uniform model a u — —2.653 (this value can be 
obtained from the value of v, quoted above, and 
the relation a u = 2 — dv). 

We have also studied the classical isotropic Heisen- 
berg model at low temperatures. Only in this limit this 
model is closed upon application of the renormalization- 
group transformation. Our procedure is appropriate to 
the study of two-dimensional systems, since only in these 
cases the non-trivial fixed point is at zero temperature. 
The renormalized parameters are [? ] : 



1 



K' A 



f m 


n 2 N 




1 ( n 3 


n 4 \ 


\K A 






' b \K A 


K B ) 



the substitution rules: 



AA---A B---B; B — > AA ■ ■ ■ A B ■ ■ ■ B . (11) 

Note that, since the model is classical, the order of the 
interactions in Eqs. (|lip is not relevant and the proporti- 
nality between the matrix of the LRGE and Ai T holds 
in this case. It is easy to show that, for K* A = K* B = oo, 
the matrix of the LRGE is given by: 



bnx/(m 
bn 3 /(n 3 



n 2 ) 2 bn 2 /(ni 
n 4 ) 2 bn±j(n 3 



n 2 f 



(12) 



Since n\ + n 2 = n 3 + n 4 = b and n 4 = in the sequences 
studied here, the eigenvalues are A 4 = n 3 /b = 1 and 
A 2 = — n 2 /b, with |A 2 | < 1. So, the eigenvalues of the 
substutitional matrix are Ai = b and A 2 = — n 2 , and the 
wandering exponnet is given by: 



lnn 2 
In b 



< 1. 



(13) 



These equations assume an aperiodic sequence built by 



Note that the value of Ai implies that v = oo and, since 
uj < 1, any aperiodic sequence with n\ +n 2 = n 3 +7ii = b 
will be irrelevant for the classical isotropic Heisenberg 
model in two dimensions. This result is supported by 
the formalims we apply in this work, since |A 2 | < 1. 



V. SUMMARY 

Witihin a real-space renormalization-group framework, 
we have studied the quantum anisotropic Heisenbeg 
model with interactions following three different aperi- 
odic sequences, on four different hierarchical lattices. We 
obtain the exact result T c = for two-dimensional lat- 
tices, while our evaluation of T c is always finite when 
d = 3. In accordance with symmetry arguments, the 
isotropic Heisenberg-model subsapce in an invariant one 
and the flow of the anisotropic models is always towards 
the Ising subspace. Our procedure allows for the calcu- 
lation of the stability of the homogeneous fixed points, 
which agrees with the Harris-Luck criterion in all stud- 
ied cases. For the relevant sequence, we established the 
presence of a new stable fixed cycle of period two and 
calculated its sepcific heat critical exponent. We also 
applied a low-temperature renormalization-group calcu- 
lation to the isotropic classical Heisenberg model in two- 
dimensional lattices: the results we obtain are exact on 
the respective hierarchical lattices and on this range of 
temperatures. We obtain that, on two-dimensional hier- 
archical lattices, all aperiodic sequences are irrelevant, in 
agreement with the Harris-Luck criterion. 
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